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Abstract 
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1 Introduction 



In the past few years, the Effective Field Theory approach has proven a useful tool to 
describe the nuclear few body problem. Like the related low energy theories of QCD, Chiral 
Perturbation Theory f^, ^ and Heavy Baryon Chiral Perturbation Theory 0, |^ in the zero 
and one nucleon sector respectively, it uses the fact that QCD exhibits a separation of scales: 
The scales associated with confinement and chiral symmetry are so high that it is natural to 
formulate QCD at low energies typical for nuclear processes only in terms of the degrees of 
freedom which emerge after quarks and gluons are confined in bound states. Below the pion 
cut , even pions are absent and nucleons emerge as the only effective degrees of freedom. A 
low energy expansion of QCD is then an expansion in powers of the typical momentum of the 
process Q in units of the high energy scale A at which the theory breaks down. This theory 
both with and without explicit pions was put to extensive tests in the two-nucleon sector 
in recent years 0-[|T^, starting with the pioneering work of Weinberg who suggested the 



usefulness of EFT's in nuclear physics [Q. It allowed for the first time for predictions which 



are systematic, rigorous and model-independent (i.e. independent of assumptions about the 
underlying non-perturbative QCD dynamics). In most cases investigated, the experimental 
agreement is within the estimated theoretical uncertainties, and in some cases, previously 
unknown counterterms describing short distance physics could be determined. Although 
in general process dependent, the expansion parameter is found to be of the order of | in 
most applications, so that NLO calculations can be expected to be accurate to about 10%, 
and NNLO calculations to about 4%. The theory with pions integrated out was recently 
pushed to very high orders in the two-nucleon sector |T3 where accuracies of the order of 
1% were obtained. It can be viewed as a systematisation of the Effective Range Theory 
(ERT) with the inclusion of relativistic and short distance effects not found in traditional 
ERT treatments. 

The simplest problem in the three-nucleon sector is nd scattering in the spin-quartet 
S wave: The absence of Coulomb interactions ensures that only properties of the strong 
interactions are probed, and the Pauli principle guarantees that the three nucleons cannot 
occupy the same point in space, making the effect of three-body forces very small. The 
phase shifts in this channel were shown to be well reproduced with and without pions both 
below and above the deuteron breakup point up to momenta of 400 MeV in the centre of 
mass frame |T^|. The explicit inclusion of pions was shown not to improve the accuracy 



even at those relatively high momenta. In the S wave, spin-doublet channel (the triton 
channel) the situation is more complicated. An unusual renormalisation of the three-body 
force makes it large and as important as the leading two-body forces. Still, by using the 
triton binding energy to fix the value of the three-body force, a reasonable description of 
the phase shift in the doublet channel can be obtained ]T7|. In this article, we apply the 



EFT with pions integrated out to the higher partial waves in the nd system, above and 
below the deuteron breakup point. All parameters are determined from A^A^ scattering, 
and three-body forces do not enter at the order we calculate, allowing one to determine the 
range of validity of the pion-less EFT without a detailed analysis of the fitting procedure. 
Since some problems depending critically on some of these channels remain unsolved, like 
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the Ay puzzle, the calculations presented are a testing ground for the applicability of the 
EFT approach to these questions. 

In Sect. I], we present a sketch of the EFT in the many nucleon sector with pions integrated 
out. The application to nd scattering in Sect. |^ is followed by a discussion of our results in 
the last Section. 



2 Formalism 

2.1 Lagrangean and Power Counting 



We now sketch the theory underlying our calculations, with more details to be found in ||16 
Three main ingredients enter in the formulation of an EFT: the Lagrangean, the power 
counting and a regularisation scheme. 

After identifying the relevant degrees of freedom in a low energy theory of QCD as the 
nucleons, one writes down the most general Lagrangean compliant with the symmetries of 
QCD. Since a theory without explicit pions is expected to break down approximately at 
the pion cut, A 100 MeV, the nucleons can be treated as non-relativistic particles with 
Lorentz invariance restored by higher order terms. Without pions, only contact interactions 
are allowed between nucleons. The first few terms in the most general, iso-spin invariant 
Lagrangean are therefore 

Cnn = NUido + -^--^)N - 
^ 2M 2M' 

- CotiN'^Pt'^Ny (A^Tp^^iV) + (2.1) 
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'{N^PfNY 


{N^PXd - 


+ ... , 







jfN) + H.c 
dfN) + H.c. 



+ 



where A^ = (^) is the nucleon doublet of two-component spinors and the sub-scripts d and t 
denote the ^Si and ^Sq channel of A^A^ scattering. For example, P^ and Pf" are the projectors 
onto the iso-scalar vector and scalar iso-vector channels, 

(Pi)l = ir^t . (P^l = (nr-t . (2.2) 

with a (r) the Pauli matrices acting in spin (iso-spin) space. A three-nucleon force appears 
at low orders only in the doublet S channel [0. The last of the nucleon kinetic energy 
terms in the Lagrangean restores Lorentz invariance ||15[. The coefficients Cj encode all 
short distance physics - like pion and uj exchanges, quarks and gluons, and resonances like 
the A - as the strengths of potentials built out of derivatives of 5 functions. In principle. 
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these constants could be derived by solving QCD or from models of short distance physics, 
but the most common and practical way is to determine them from experiment. 



Since the Lagrangean (|2.1| ) consists of infinitely many terms only restricted by symmetry, 
predictive power is ensured only by the second ingredient of an EFT: a power counting 
scheme, i.e. a way to determine at which order in a momentum expansion different contri- 
butions will appear, keeping only and all the terms up to a given order in calculations. The 
dimensionless, small parameter on which the expansion is based is the typical momentum 
Q of the process in units of the scale A at which the theory is expected to break down, 
e.g. in a pion-less theory in units of the pion mass. Assuming that all contributions are of 
natural size, i.e. ordered by powers of Q, the systematic power counting ensures that the 
sum of all terms left out when calculating to a certain order in Q is smaller than the last 
order retained. This way, an EFT allows for an error estimate of the accuracy of the final 
result. 

In the two-nucleon sector, finding such a power counting scheme is complicated by the 
fact that the scattering lengths of both the ^Sq and ^Si channel in A^A^ scattering are 
unnaturally large {1/at = —8.3 MeV, l/a^ = 36 MeV) because of the existence of a 
shallow virtual and real bound state, the latter being the deuteron with a binding energy 
B = 2.225 MeV and hence a typical binding momentum 7^ = yMB ~ 46 MeV. This 
lies well below the expected breakdown scale of the EFT and suggests that the system 
is close to a non-trivial fixed point^. Due to this fine tuning, the naive low momentum 
expansion has a very limited range of validity {Q < l/at,Q < l/fld). The alternative is 
to expand in powers of the small parameter Q/A ^ 1 but to keep the full dependence on 
Qott, Qdd ~ 1- This is clearly necessary if one is to be able to describe bound states within 
the effective theory. Some interactions have therefore to be treated non-perturbatively in 
order to accommodate these bound states. In the effective theory, the fine tuning present in 
the two-body interactions arises as cancellation between loop and counterterm contributions. 
This makes it tricky to do the power counting in a naive approach since it cannot be applied 
diagram by diagram, but only to whole classes of diagrams which are furthermore not 
always easily identified |jl9|. A convenient way of dealing with this problem was suggested 
in [^]. It consists in using dimensional regularisation and a new subtraction scheme (Power 
Divergence Subtraction, PDS) in which not only the poles Feynman graphs exhibit in 4 
dimensions are subtracted but also the poles in 3 dimensions. This regularisation procedure 
is chosen to explicitly preserve the systematic power counting as well as all symmetries at 
each order in every step of the calculation. If one chooses the arbitrary scale /i arising in 
the subtraction to be of the order /i ~ Q ~ l/a^ ~ l/a^, one finds 



(-1) 







MQ 



^See e.g. ^] for an analysis of this short distance fine tuning problem in terms of eritieal points of the 
renormalisation group. 
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Co 



MAQ2 



for the short distance coefficients of both the ^Si and ^Sq channeL Here, as in the following, 



we split the coefficients of the leading four point interactions, Co, into a leading (Cg and 
a sub- leading piece (C^^), where the super-script in parenthesis denotes the scaling of the 
coefficient with Q. 



Using ( |2.3| ) one can estimate the contribution of any given diagram. As will be demon- 
strated in the next sub-section, it is found that for nucleon-nucleon scattering, the leading 
order contribution is given by an infinite number of diagrams forming a bubble chain, with 
the vertices being the ones proportional to Cq. Higher order corrections are perturbative 
and given by one or more insertions of higher derivative operators. In the three-body sector 
discussed in this paper, even the leading order calculation is too complex for a fully ana- 
lytical solution. Still, the equations that need to be solved are computationally trivial and 
can furthermore be improved systematically by higher order corrections that involve only 
integrations, as opposed to many-dimensional integral equations arising in other approaches. 

Following a previous paper on quartet S wave scattering in the nd system [0, we choose 
a Lagrangean which is equivalent to but contains two auxiliary fields rf* and with 
the quantum numbers of the deuteron and of a di-baryon field in the ^Sq channel of A^A^ 
scattering respectively, such that the four nucleon interactions are removed: 



Nd 



^ ° 2M 2M^ 



;i5o + 
(i9o + 



AM' 



+ [d'\N^P',N) + H.c] + 
+ yt [t^\N^PfN) + H.c] + . 



(2.4) 



It was demonstrated in [|l6l that this Lagrangean is more convenient for the numerical 



investigations necessary in the three-body system. Analogously to the Cq's, we split the A's 
into leading (A*^"^)) and sub-leading pieces (A*^°^). Performing the GauBian integration over 
the auxiliary fields and in the path integral followed by a field re-definition in order to 
eliminate the terms containing time derivatives, we see that the two Lagrangeans (|2.1|) and 
are indeed equivalent up to higher order terms when one identifies 



A(-i) 
A(o) 



(-1) 



Cc 
MC2 

MCo ' 



(2.5) 



1)n2 



{CD 

MCo 
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for each of the two S wave channels separately. The "wrong" sign of the kinetic energy 
terms of the auxiliary fields therefore does not spoil unitarity. The new coefficients scale 
from (p. 31) and (12.51) in each channel as 



M 



A(°) ~ g , (2.6) 



A 

M2 



2.2 The Two-Nucleon System to NNLO 

It is useful to consider the solution to the two-body problem before turning to the three-body 
problem. As the ^Si and ^Sq channels do not differ in their power counting, we concentrate 
on the deuteron channel. The considerations in the ^Sq channel are analogous. More details 
about the EFT without pions in the triplet channel can be found in 

Any diagram can be estimated by scaling all momenta by a factor of Q and all energies 
by a factor of Q'^/M. The remaining integral includes no dimensions and is taken to be of 
the order and of natural size. This scaling implies the rule that nucleon propagators 
contribute one power of M/Q'^ and each loop a power of /M. Relativistic corrections 
to the kinetic energy of the nucleon scale like jj^ and hence only enter as insertions into 
the non-relativistic nucleon propagator at NNLO in Q, suppressed by additional powers 
of 1/M. The vertices provide powers of Q according to (|2.6| ), implying that the deuteron 
kinetic energy term is sub-leading compared to the A^"^'' term. Thus, the bare deuteron 
propagator is just the constant — i/A^ ^\ Using these rules, there is an infinite number of 
diagrams contributing at leading order to the deuteron propagator, as shown in Fig. |l|, each 
one of the order 1/[MQ). The linear divergence in each of the bubble diagrams shown in 



Figure 1: The deuteron propagator at leading order from the Lagrangean (\2.^ ). The thick 
sohd hne denotes the bare propagator — the double line its dressed counterpart. 

Fig. |l] does not show in dimensional regular isat ion as a pole in 4 dimensions, but it does 
appear as a pole in 3 dimensions which we subtract following the PDS scheme. 

The full leading order propagator iA^^(p) of the deuteron field consists hence of the 
geometric series shown in Fig. |T], 



6 



where /i is the arbitrary renormahsation scale introduced by the PDS scheme, showing that 
^VyJ is renormahsation group dependent, too. Physical observables like amplitudes are 
of course independent of the choice of /x, as will be demonstrated now. We determine the 
free parameters by demanding that the deuteron propagator has its pole at the physical 
deuteron pole position as 

Since 7^ ~ Q, the choice fi ^ Q makes the coefficient CqJ^^ indeed scale as required in the 
power counting scheme ( p.3|) . The deuteron and di-baryon propagators are hence in terms 
of physical quantities at leading order given by 



iA-(,) ^ A.(.) ^ ^ . (2.9) 

It - y^- Mpo - le 

The typical momentum 7^ = l/a^ of the virtual bound state is extracted from the scattering 
length in this channel. Each deuteron or di-baryon propagator is accompanied by a power 
ofl/{My^Q) = M/{AQ). 

The kinetic energy insertion and the mass insertion proportional to A^^"* are suppressed 
by one power of Q according to the scaling properties (|2.6| ) and hence enter at NLO. We 
choose to keep the A^^^^'s to be the same as at LO, so that (p.8|) is still valid at NLO 



and higher order calculations will not necessitate a re-fitting of lower order coefficients ^1 



The remaining pieces, are then parametrically smaller (as indicated by ( |2.6D ) and are 
included perturbatively. Two conditions per partial wave are now necessary to fix the new 
constants at NLO. One condition is that the pole position of the bound state does not 
change. Fig. S 

2 

Ai°) = ^ (2.10) 

The second condition in the "^Si channel can be that the deuteron pole has the correct 
residue Zd ||22| , 

^ = ^.-1 , Z^' = l-^dP, , (2.11) 
where pd is defined by the effective range expansion around the deuteron pole, 

k cot 6 = --fd + ^Pd{k^ + Id) + ■■■ , (2.12) 

for the nucleon-nucleon phase shifts. In the singlet channel, no real bound state exists 
whose properties have to be described, so we impose the condition that the effective range 
expansion 

kcot6 = -— + -rotk^ + ... (2.13) 
at 2 
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.(-1) 



Ai°^ and yd, l^/^l^). The dressed 



Figure 2: The first condition on the coefficients A 
deuteron propagator has its pole at the physical binding energy at LO, and higher order 
corrections do not change its position. The cross denotes a kinetic energy insertion, the dot 
an insertion of A^?^ . 



is reproduced. Solving (|2lO| / |2lIl / ^l3D with (|2j) leads to 



_2_ IJ'-ld 
M 



A 



(-1) 



Pd 



2 12 

M ~ 



It 



rot 



A 



(0) 



7i ^ 
M ' 



(0) 







(2.14) 



Vd 



Id 



M2 Zd 



Vt 



Stt 1 

(-2) 



In a NNLO calculation, C2 is split like Co, with C2 still given by the expression above 



C2 a NNLO correction parametrically small against C2 ^' ■ Likewise, Co has now a NNLO 
component Cq^'* which is again fixed by the condition that the deuteron or di-baryon pole 
position is unchanged. This way, higher order calculations do not necessitate a re-fitting of 



(-2) 



lower order coefficients |2^. In addition, four new interactions enter at NNLO [|T5[: The 
first is a term which induces SD mixing in the two-nucleon system and will be shown to be 
irrelevant for our problem. Relativistic corrections to the nucleon propagator, — 2M ~ 
enter at NNLO in the Q counting, but are suppressed by two more powers of M and hence 
can be neglected. Finally, there are A^A^ interactions containing four derivatives. In the 
triplet channel, for example, a term 

- 4 {N^Plid - dfNy {N^P'S - dfN) (2.15) 
lb 

enters in the Lagrangean (|2.1|) which by GauBian integration and field re-definition can be 
shown to be equivalent to a term in the modified Lagrangean (|2.4|) of the form 



Ma 



id 



c. 



(-2) 



ido + 



AM 



n 2 



(2.16) 



However, it was shown in [|T^ that the renormalisation group equations determine the 
coefficient C4 through the lower order constants Cq and C2, i.e. that the value of this 
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constant is determined already from the NLO parameters as [22| 



G 



4d 



(2.17) 



This is a generic feature, the leading piece of these higher derivative terms are determined 
by the LO and NLO coefficients. Their independent contribution start at N3L0 but can be 
re-summed into a full deuteron (di-baryon) propagator 



iA^^ (p; N3L0) 



iAr(p;N3L0) 



i5'^ Ad{p; N3L0) 
47ri 



-Id . 2 

i6^^Atip; N3L0) 



ip,(Mpo - f + ij) + a/ ¥ - Mpo - is 



My, 



-It 



lr,tiMpo-^) + J^-Mpo-ie 



(2.18) 



3 The nd System 
3.1 Quartet Channel 

Turning to the nd system, we first consider the quartet channel as its treatment is analogous 
to the doublet channel without the difficulties of a coupled channel equation. For details, 
we refer again to the treatment of the quartet S wave in ||16|| . 

All three nucleon spins are aligned, so that only properties of the triplet A^A^ system 
are probed. It has been shown in [|T^ that all diagrams in which three nucleons interact 
only via the LO deuteron-nucleon potential yd are of the same, leading order K/{MQ'^). 
Summing all "bubble-chain" sub-graphs of Fig. |l| into the deuteron propagator, one arrives 
at the integral equation pictorially represented in Fig. 0. To derive this equation for the 





Figure 3: The Faddeev equation (jSTlp which needs to be solved for at LO in the quartet 
channel. 



half off-shell amplitude at LO, let us choose the kinematics as in Fig. ^ in such a way that 
the incoming (outgoing) deuteron line carries momentum k (p), energy A;^/(4M) — 7J/M 
(fc^/ (4M) — + e) and vector index j (i). The incoming (outgoing) nucleon line carries 
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momentum —k {—p), energy fc^/ (2M) (fc^/ (2M) — e) and spinor and iso-spinor indices a and 
a {(3 and h). Hence, the incoming deuteron and nucleon are on shell, and e denotes how far 

the outgoing particles are off shell. We denote by (itoq)|^^ (k, p, e) the sum of those diagrams 
with the kinematics above forming the half off-shell amplitude in the quartet channel, and 
read off from the lower line of Fig. ^ the integral equation: 

(iti^y {k,p,e) = yl{a^a'f,6',—^ , = + 

4M M ^ 2M '■^ 

+ f (-'-r, ii I ^ (oi: Ck, ?. < + *) X (3.1) 

2M ^ ^0 2M ^ ^ 



X 



'4M ~ M + ^^ + 10 2W- + 1^ 



The integration over go picks the pole at go = (^^ ~ ^^)/(2M) — e + ie. After that, we 
set e = (/c^ — p^)/(2M), integrate over the angle between A; and p weighted with the 

Legendre polynomial of the first kind Pi{k ■ p) to project onto the Ith partial wave and set 
i = {1 + i2)/v^, j = (1 — i2)/\/2, a = /3 = 1, a = b = 1, to pick up the spin quartet part. 
Denoting this projected amplitude by itoq(/c,p), we find 

oo 



where E = 3k'^/{AM) --fl/M is the total energy and the Legendre polynomials of the second 
kind Qi are defined as 

1 

Q,{a) = \ ! dx ^ . (3.3) 
2 1 X + a 



This Faddeev equation was first derived in p3[. Notice that when p = k, all external 



legs are on-shell. Although the pole on the real axis due to the deuteron propagator is 
regulated by the ie prescription and the logarithmic singularities occurring above thresh- 
old are integrable for each partial wave, both cause numerical instabilities. We used the 



Hetherington-Schick p4| , p5| , [26[] method to numerically solve ( |3.2| ). The basic idea is to 
perform a rotation of the variable q into the complex plane by an angle large enough in order 
to avoid the singularities in and near the real axis but small enough not to cross the singu- 
larities of the kernel or of the solution. One can show that the singularities of the solution 



are not closer to the real axis than those of the inhomogeneous. Born term [27]. The solution 
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on the real axis can then be obtained from the solution on the deformed contour by using 
( |3.2| ) again, now with k and p on the real axis and q on the contour. The computational 
effort becomes trivial, and a code runs within seconds on a personal computer. 

To obtain the neutron-deuteron scattering amplitude, we have to multiply the on-shell 
amplitude toq(fc,p) by the wave function renormalisation constant. 



(0, 



Oq(^) 
1 



I Zq t, 

d 



1 



MSI 



P-- 



(3.4) 



In contradistinction to itQq(/i;,p), the scattering amplitude iTQ^{k) depends on the parame- 



ters i/d and Ai only through the observable 7^. 



At NLO, we have additional contributions only from modifications of the deuteron, as 
Sect. p.2| demonstrated, namely the deuteron kinetic energy and A^°^ insertions depicted 
in Fig. 0. One might note that the diagrams are all finite thanks to the re-formulation of 




Figure 4: The NLO contributions to nd scattering in the quartet channel. Notation as in 
Fig. |. 



the original Lagrangean into (p.4| ) and hence can be treated numerically fl^. Each 
contribution is inserted only once. Because of the absence of interactions which mix different 
partial waves, the NLO contributions are given by 



2M ^ '■^ 



where the corrections to the deuteron are 



The go integration picks the nucleon pole resulting in a one dimensional integral 



(3.5) 



(3.6) 



27r2 



iA„ 



q- 



-\ 2 



^4M M 2M' 



X 



(3.7) 
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X iXv 



'AM M 2M' 



which has to be performed numerically since t'l^l^{k,p) is known only numerically. 
Finally, the wave function renormalisation constant Z at NLO is found from 

1 



1 
Z 



Zq + Z\ 
. d 



Zq Zq 

1 



dpo iArf(p) + iAd{p)iId{p)iAd{p) 



-7^ 

p=0 



Zo opo 



'd 
p=0 



(3. 



as 



Z, = Z^ . (3.9) 
The NLO amplitude of the Ith partial wave in the quartet channel is therefore given by 

Tj;\k) = zt(i\k,k) 



{Zo + z,){S{k,k) + mk,k)) 



(0, 



T^lik) + ZotZik. k) + z,tZ{k. k) = n^'ik) + rrm . 



,(0 



(0/ 



(3.10) 



The phase shift for each partial wave is extracted by expanding both sides of the relation 



rf(A;)^To«(A:) + r,(')(fc) = :^ 



(0/ 



37r 



1 



M fccot5W(fc) - ik 
in Q with d^''^ = 5^^^ '■'^ + 6^^^ *■'•* + ... and keeping only linear terms, 

2i V 

kM T^''{k) 



(3.11) 



1 , / 2iA;M 
- n 1 - 
2i V 



ovr 



3n (l + ?^T^l^k) 



(3.12) 



At NNLO, there are two kinds of contributions. One comes from the same deuteron kinetic 
energy and A^^^ terms appearing at NLO, but now inserted twice. The other comes from 
inserting an operator describing the tensor force which appears for example in the spin 
triplet channel of the Lagrangean (|2.1|) as 



( Pi(9 - d)\d - dV - ^P^id -dr]N 



+ H.c. (3.13) 



and is generated by the term 

'-Nd — hid 



N[{d- d)\d - dY - -s'Kd - df ] P^N 



H.c. 



(3.14) 
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in the Lagrangean (p.4|) after the GauBian integration over the deuteron field is performed. 
This operator mixes spin and angular momentum and hence produces a splitting and mixing 
oft amplitudes with the same spin S and angular momentum L but different values of the 
total angular momentum J. This splitting is very important in determining spin observables. 
However, although our calculation for the phase shifts themselves is at NNLO and should be 
precise, the splittings start only at this order and should not be very realistic. We postpone 
the analysis of these splittings to a future publication where the spin observables will also be 
addressed and limit ourselves here to the phase shifts averaged over different values of J. At 
NNLO the tensor force appears only linearly, and a somewhat lengthy but straightforward 
calculation shows that, after this average, its contribution vanishes. 

The following trick allows for a very simple calculation of the double insertion of the 
kinetic and A^^^ operators. If they were included to all orders, the only change from the LO 
calculation would be the use of the NNLO deuteron propagator ( p.l8| ) in instead of 
the LO one, (|2.9| ). This procedure re-sums some contributions of order N3L0 and higher, 
but not all of them, so it can be justified only up to NNLO. Therefore, a spin-averaged 
NNLO calculation is simply performed by replacing the LO deuteron propagator (|2.9|) by 
the propagator ( |2.18| ) which summarises all effective range effects, in the integral equation 
(|3.2|) . As the resulting amplitude is unitary, is used directly to extract the phase shift. 

This procedure has already been used in the calculation of the S wave quartet scattering 
in ill. 



3.2 Doublet Channel 

The calculation for the spin 1/2 channels proceeds in an analogous way to the one in the 
spin 3/2 channels discussed above. The main difference is that the spin zero di-baryon field 
t also contributes in intermediate states of nd amplitudes. The analogue of ( p.l|) is then 
a system of two coupled integral equations for the Nd Nd amplitude it^^ Nd^Ndi^^P^ ^) 
and for the Nd — > Nt amplitude it}^ jsid-^Nti^^V^ pictorially represented in Fig. ^: 





+ 




Figure 5: The coupled set of Faddeev equations ( |3. J5|/|3n^ ) which need to be solved for 
tod at LO in the doublet channel. The double dashed line denotes the di-baryon field t"^. 
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4Af A/ ^ 



2AI 



+ 



(2^ 



^2 2 
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After projecting on the spin 1/2 channel using 
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and following the same steps as in the quartet case one finds 
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The modifications for tlie NLO and NNLO calculation are straightforward. Fig. ^ shows 
the graphs giving the NLO corrections. At NNLO, the propagators iAd/t{p) in the above 
integral equations are replaced by iAii/t{p; N3L0) of ( |2.18| ). 







-iA_ 

:::0:: 



(0) 




Figure 6: The NLO contributions to nd scattering in the doublet channel. 



4 Discussion and Conclusions 

We used he = 197.327 MeVfm, a nucleon mass of M = 938.918 MeV, for the NN triplet 
channel a deuteron binding energy (momentum) of -B = 2.225 MeV (7^ = 45.7066 MeV), 
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a residue of = 1.690(3), and for the A^A^ singlet channel an ^Sq scattering length of 
at = —23.714 fm and effective range vqi = 2.73 fm. In Fig. we present the results of our 
computations for the real and imaginary parts of the quartet partial waves / = to 4 in the 
centre-of-mass frame. Fig. ^ shows the real and imaginary parts of the doublet partial waves 
/ = 1 to 4. The ^Sa partial wave was already computed in a previous publication [|l^ with 
a choice of parameters and method to extract the phase shift which differs from the one 
used here formally only at higher order in the power counting, and indeed the difference is 
marginal. Comparison of the LO with the NLO and NNLO result demonstrates convergence 
of the effective field theory, with the expansion parameter found to be about |. We therefore 
claim that the NNLO calculation has an error bar of less than 4 percent. It is interesting 
to note that the NLO correction is for high enough energies sometimes sizeable, while the 
NNLO correction is in general very small. In contradistinction to the NN case, the hierarchy 
the phase shifts exhibit at large momenta with increasing / is not due to suppression of the 
higher partial waves in powers of Q, as the diagrams scale for all partial waves in the same 
way. Rather, it is easy to see that especially for the quartet channel, already the behaviour 
of the LO integral equations ( |3.2| ) and ( |3.18| / |3TT9|) as p = /c S> ■jd/t makes the Ith partial 
wave asymptotically go like the /th Legendre polynomial of the second kind Qi{j), and 
Qiiz)/Qi-iij) ~ 0-33 — 0.45, which is by accident close to the expansion parameter Q. 

Below the deuteron breakup point, our calculation is in good agreement with the vari- 
ational calculation of Kievsky et al. which is based on the AV18 potential. We also 
compare to partial wave analyses of pd scattering experiments by Huttel et al. |^ and 
Schmelzbach et al. as Coulomb and chiral symmetry breaking effects can be neglected 



at high momenta. Above Ecm ~ 15 MeV, i.e. a cm-momentum of about 140 MeV, one 
expects the pion-less theory to diverge from experiment because not later than then should 
the pion manifest itself as an explicit degree of freedom of the low energy theory. 

Within the range of validity of this pion-less theory we seem to have good convergence 
and our results agree with potential model calculations within the theoretical uncertainty 
of (1/3)^ ~ 4%. That makes us optimistic about carrying out higher order calculations of 
problematic spin observables where our approach will differ from potential model calcula- 
tions due to the inclusion of three-body forces. For example, the nucleon-deuteron vector 
analysing power Ay in elastic Nd scattering at energies below E'lab ~ 30 MeV does not have 
a satisfactory description in any realistic potential model of A^A^ scattering, see [^, |3^ for 
details. It remains to be seen whether the EFT approach can be of help there. 

We close with a note comparing the EFT approach at this order to realistic potential 
model calculations. Since we use as input parameters to NNLO only the scattering lengths 
l/7(i, cit and effective ranges pd, rot of the two body system in the S wave channels, we ex- 
pect any realistic potential model which reproduces these numbers to agree with our results 
within the error bar of an NNLO calculation, ~ 4 %. The absence of counterterms at this 
order (e.g. of three body forces) guarantees that an effective theory of a realistic potential 
model will be identical to the EFT discussed here. Hence, a potential model can already 
serve for the purpose of comparing to our results at energies considerably higher than the 
deuteron breakup point. However, the EFT approach is systematic and rigorous, allowing 
for error estimates. The small number and simplicity of the interactions in the Lagrangean 
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at NNLO makes the computations considerably simpler: The LO integral equation is only 
one-dimensional and can be obtained in very short computing time, and at NLO, only (par- 
tially analytical, partially simple numerical) one dimensional integrations are necessary, in 
contradistinction to many dimensional integral equations in traditional approaches. Calcu- 
lations of phase shifts in any approach, be it potential models or EFT with pions as explicit 
degrees of freedom, are hence bound to reproduce our results within error-bars in the range 
of validity of the pion-less theory, allowing for valuable cross-checks. Finally, rigorously 
determined three body forces will enter at higher orders of the EFT approach and lift the 
similarity to potential model calculations. 
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Figure 7: Real and imaginary parts of the Grst five partial waves in the quartet channel of 
nd scattering versus the centre-of-mass energy. The dashed line is the LO, the dot-dashed 
the NLO and the solid line the NNLO result. The calculations ofKievsky et al. (crosses) are 
from Refs. below breakup (Ecm = B) and above breakup. The phase shift analyses 
of pd data by Huttel et al. / |30|/ and Schmelzbach et al. J3l]/ are presented as open squares 
and diamonds, respectively. 
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G Wave 
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Figure 8: Real and imaginary parts of the first four liigher partial waves in the doublet 
channel of nd scattering versus the centre-of-mass energy. Notation as in Fig. ^ 
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